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Abstract. We propose a new way to detect individ- 
ual briglit Ultra-Higli Energy Cosmic Ray (UHECR) 
sources above background if the Galactic Magnetic 
Field (GMF) gives the main contribution to UHECR 
deflections |1|. This method can be directly appUed 
to maps given by experiments. It consists in start- 
ing from at least two high energy events above 
6 X 10^^ eV, and looking at lower energy tails. We 
test the efficiency of the method and investigate 
its dependence on different parameters. In case of 
detection, the source position and the local GMF 
deflection power are reconstructed. Both reconstruc- 
tions are strongly affected by the turbulent GMF. 
With the parameters adopted in this study, for 68% 
of reconstructed sources, the angular position is less 
than one degree from the real one. For typical 
turbulent field strengths of 4 fiG at the Earth position 
and 1.5 kpc extension in the halo, one can reconstruct 
the deflection power with 25% precision in 68% of 
cases. 

Keywords: ultra-high energy cosmic rays, magnetic 
fields. Galaxy. 

I. Introduction 

The HiRes experiment first observed the cutoff in the 
spectrum fT\. The UHECR sources are still unknown. 
Some possible candidates can be studied from the point 
of view of acceleration mechanisms -see for example 
Refs. f3l, m, fSl. 

Their detection is difficult due to the UHECR deflections 
in the extragalactic and Galactic magnetic fields. One 
can either wait for enough data at the highest energies 
to find a class of UHECR sources, or search for the 
first individual brightest sources. In this study, we 
investigate this second possibility. 
We propose here a method to find Ultra-High Energy 
(UHE) proton or light nuclei sources on top of 
background and reconstruct their positions. We show 
that one has to start from at least two events with 
energies above 10^^'^eV^ 6 x 10^^ eV before looking 
at lower energy tails, in order to avoid confusion by 
the background. We investigate the performance of this 
method depending on the different relevant parameters. 
We assume as in [6| that UHECR deflections due 
to the extragalactic magnetic fields are negligible 
compared to those due to the GMF. For an example 
of non negligible extragalactic contributions, see |7|. 
The GMF is divided in two components: regular and 
turbulent. For the regular one, we take the Prouza 
and Smida model (SI, ||9l. It consists of a thick disk. 
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Fig. 1 : Projections of the arrival directions of cosmic rays (CR) 
emitted by a source (cross) in the plane tangent to the celestial 
sphere, and centered on the highest energy event (filled square 
in (0,0)) emitted by the source. Deflections X and Y measured 
on two orthogonal axes, in degrees. Open symbols for source 
CR deflected by the regular GMF only. Filled symbols for 
the same CR when the turbulent component is also taken into 
account. Shapes correspond to CR energies (triangles: between 
10^^'' and lO^'''^ eV, circles: between lO^^'^ and lO"'' eV, 
and squares: above 10^^'' eV). 



and of toroidal and dipolar fields. Following the latest 
knowledge on the GMF ifTOl . ifTTI . we modify some 
parameters and take an exponentially decaying profile 
along the Galactocentric radius. The field strength close 
to the Sun is Bj-^g = 1 — 3/LtG (currently admitted 
strength: Breg — 2 /iG). The model for the turbulent 
component is described in [T|. The root mean square 
(RMS) value of the turbulent field decays as 1/r along 
the Galactocentric radius r, and exponentially along the 
direction orthogonal to the Galactic plane. The plane 
thickness is set here to 1.5 kpc. The turbulent field 
strength close to the Sun is Sturb = 2 — 8 /iG (currently 
admitted strength: i?turb — 4 fiG). 

For source events, we assume power law acceleration 
spectra, with a power law index equal to 2.2 and a 
maximum energy set to 10^^ eV. In most cases, we take 
a proton source and set its distance to the Earth to 
50 — 100 Mpc. Protons are propagated to the observer by 
using the results of Ref. |12| for energy losses, which 
creates a "bump" in the source spectrum ||T31 . Then, the 
generated cosmic rays (CR) are deflected in the regular 
GMF (see e.g. 1 14 1), and in the turbulent GMF. For more 
details on turbulent GMF deflections, see Refs. ifTSl . fl\. 
The effects of both components are shown in Fig. [T] 
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The regular component deflects the CR along a curve 
and the turbulent one spreads them "randomly" around 
it in a sector-shaped region. When there are many events 
with the same energies, they are not spread uniformly, 
as discussed in References 1141 . 1151 . [16J. 
Results depend on many parameters of the model and 
on the considered location on the sky. However, in 
practice, there are only two essential parameters: The 
local deflection powers -see Ref. [[T]- of the regular GMF 
(denoted D) and of the turbulent GMF. 
Background events are simulated according to the energy 
spectrum measured by HiRes IS] and the exposure of 
Telescope Array \Yl\ . 

II. Method 

The following work is done on sky maps which 
only display events with energies bigger than a given 
threshold, ii'th- The idea is to look for tails of lower 
energy events around at least two nearby events with 
energies E > 10^^-^ eV. 

We start by taking a circle around an event with 
El > 10^^'^eV. Its radius R can be optimized and 
the value we take depends, among others, on Eth 
and D. We will call R an "internal" parameter of the 
method. The tail of lower energy events is searched by 
assuming that all events from the source are located in 
a sub-region of the circle, that has a sector shape with 
a central axis given by the second highest energy event. 
The energy E2 < Ei of the second highest energy event 
is above a given threshold E2 > i?2min, and its distance 
to the highest energy event must be compatible with an 
emission from the same source. This distance should be 
lower than PD/E2 — D/Ei, where /3 is another internal 
parameter, and D is the value initially assumed for the 
local deflection power. Fig. |2] shows the probability to 
detect the source and the probability to be confused 
by the background, for different values of £'2min- One 
can see that one has to require £'2min = 10^^'^ eV, in 
order not to detect some background more frequently 
than the source. That is why we search for at least two 
nearby events at the highest energies, E > 10^^'^ eV. 

Then, as a discriminator of source events, we use 
a correlation coefficient -see Ref. |16|. One takes the 
energies E and the coordinates X' of the events lo- 
cated in the sector, where X' is the axis containing 
the first event and the center of mass of the other 
considered events. Each event is tested and if removing 
it increases the correlation coefficient C{1/E,X'), it is 
then definitely removed. Otherwise, it is definitely kept. 
If C{l/E, X') is finally larger than a given value, Cmin, 
there is detection: Detection of the source in case of an 
initial source doublet or confusion by the background in 
case of an initial background doublet. 
Ref. im discusses the optimization of the internal pa- 
rameters. 

In case of source detection, we finally reconstruct its 
position along the X' axis, by taking the remaining 
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Fig. 2: Dependence on £'2min, the minimum energy for the 
second highest energy event. Thick solid line for the probabil- 
ity to detect the source. Thin solid line for the probability 
to detect some background. Results presented for a source 
of luminosity above 10^^'* eV equal to 2.6% of the total 
luminosity in UHECR, Eth = 10^^'^ eV, 5000 events above 
10^'^ eV in the whole sky, B^cg = 2/iG and Bturb = 4pG. 



events and fitting 1/E versus X' with a straight line. 
This also gives an estimate of the local regular GMF 
deflection power, D. 

III. Study of the parameter space 

We study below how the efficiency of the method 
presented above changes with values of physical or 
experimental parameters. 

There are two requirements for source detection. First, 
the probability to detect the source doublets has to be 
large enough, in order not to remove too much source 
signal. Second, the probability to be confused by some 
background has to be low enough compared to the 
probability to detect the source signal. 
We define here a "cluster" as a group of at least three 
nearby events with very high energies, E > 10^^'^ eV. 
The probability to have a cluster of background events 
is very low in case of reasonable statistics and no 
-or low- anisotropy. For example it is ~ 1%, with 
5000 events in the whole sky above 10^^ eV and no 
anisotropy. On the contrary, the probability to have a 
background doublet of very high energy events is, in 
most cases, not negligible compared to the probability 
to have a source doublet. That is why we separate 
the two cases -doublet or cluster- below. In case of a 
doublet, we apply the method depicted in the previous 
section, whereas in case of a cluster, we do not need it 
since the source is already detected in ^ 99% of cases. 
Applying the method reduces the background more 
than the signal from the source, but its side-effect is 
to reduce the signal from the source in any case. That 
is why it is not worth applying it when the probability 
to be confused by some background is already very 
low. However, when statistics increase, or in case of a 
large anisotropy above 6 x 10^^ eV, background clusters 
become more frequent and the method would then also 
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Fig. 3: Dependence on the regular GMF strength, Breg, with a 
constant ratio -Brcg/Stmb set to 0.5. Eth = 10^^'^ eV. Source 
luminosity and number of events on the sky as in Fig. [J] Dotted 
line for the probability to have a doublet or a cluster from the 
source. Dashed-dotted and dashed lines for the probabilities to 
have a source or background doublet, respectively. Same key 
for solid lines as in Fig. [J] 
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Fig. 4: Dependence on the turbulent GMF strength, Bturb, with 
a constant regular field set to Bicg = 2pG. Same values for 
the other parameters and same line types as in Fig. |3] 



be interesting for clusters of 3 events. 

The figures |3] and |4] show the dependence of results 
on the regular and turbulent GMF strengths. Fig. [3] 
is computed for different strengths of B^cg between 1 
and 3 /iG, and with a constant ratio between regular 
and turbulent field strengths. For Fig. H) B^cg = 2 /iG 
and the turbulent field strength i?tuib varies between 2 
and 8 fiG. The internal parameters of the method have 
been optimized for each field strength. The dotted line 
is the probability that a proton source of luminosity 
above 10^^^ eV equal to 2.6% of the total luminosity 
in UHECR emits a doublet or a cluster of events with 
energies E > lO^^'* eV. The dashed-dotted and the 
dashed lines are respectively the probabilities to have a 
source or background doublet. The thick and the thin 
soUd lines represent the probabilities to respectively 
detect the source or the background. One can notice 



that even for the largest field strengths in these ranges, 
the ability of the method to detect the source is still 
quite acceptable. The difference between the dashed and 
the thin solid lines shows the efficiency to remove the 
background. The difference between the dashed-dotted 
and the thick solid lines is smaller: The source signal 
is less removed than the background. The difference 
between dotted and dashed-dotted lines corresponds, in 
practice, to cases when the source is already detected 
thanks to clusters of events. 

These results are computed for an isotropic sky, even 
at the highest energies. We checked that for reasonable 
values of anisotropies (if all events with energies above 
JQ19.8 gY ^j-g located in a fraction of the sky larger than 
25%), results would be affected in a linear way. For 
example, in the case of a 25% fraction, the probabilities 
to have a background doublet and to detect it would be 
both approximately multiplied by 4. 
We have studied the dependence on source parame- 
ters. We show in Ref. yj results for the dependence 
on the source luminosity. We point out that with a 
regular GMF deflection power, for example close to 
_D ~ 5° X 10^^'^ eV for protons, one can still detect 
a light nuclei source, but only by looking for clusters if 
the source luminosity is sufficient. In this case, the ratio 
"signal/background" for doublets would be too small to 
detect such sources. We have also presented results on 
the dependence on experimental statistics. There is no 
clear and no general best energy threshold Eth to look 
for lower energy tails of events. When Eth is decreased, 
the probability to detect a source doublet through the 
method increases, but the probability to be confused by 
a background doublet increases with a larger rate. For 
the example considered in Ref. [Tl, Eth — 10^^'^eV~ 
4 X 10^^ eV is a good compromise between sufficient 
source detection and sufficient background rejection. 
The impact of the total amount of experimental data 
mostly affects the way the detection should be done, 
depending on the source luminosity. For example, with 
a bright source and large statistics like 10'* events above 
lO^^eV, clusters of source events appear more often 
and looking for them becomes a more efficient way to 
track the source than applying the method for doublets. 
However, when the amount of data still increases, the 
method used previously for doublets will also start to 
be needed for clusters, so as to reject background ones. 

IV. Reconstruction of the source position 
AND of the regular GMF deflection power 

In this section, we discuss the ability of the method 
to reconstruct the source position and the local regular 
GMF deflection power D. 

With the assumptions considered here, and Bturb = 
4/uG, the position of a source located in a region of 
the sky where D ^ 5° x 10^^'^ eV can be reconstructed 
with an accuracy of 1 degree in about 68% of cases - 
see distribution with thick solid line in Fig. |5] Thus, the 
angular resolution of this method for the source position 
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Fig. 5: Distribution of distances, in degrees, between the 
reconstructed sources and the real source. The source is located 
in a region of the sky where D ~ 5° x 10^®'® eV. Thin solid 
line for no turbulent component (Bturb ~ 0/iG), dashed line 
for Bturb = 2 fiG, and thick solid line for Btmb ~ 4 pG. 
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Fig. 6: Distribution of reconstructed deflection powers in ° x 
]^gi9.6 gY Same parameters and line types as in Fig.jS] 



reconstruction is of the order of the experimental angular 
resolution. 

Another noticeable result is that, in the same conditions, 
the local deflection power D of the regular GMF com- 
ponent can be reconstructed with up to 25% precision in 
68% of cases, as shown in Fig. |6] -see thick solid line. 
The above results are mainly spoiled by the turbulent 
field. For comparison, we put in the same plots the 
distributions with i?turb = O/iG (thin solid line) and 
Sturb = 2 /iG (dashed line). The lower the turbulent 
field strength, the better the precisions on both results. 
For these lower Bturb strengths, the distance on the sky 
between the reconstructed and the real sources would be 
mostly affected by the experimental angular resolution 
and not by the method. 

The strength of the regular field (for i?rcg = 1 — 3 /iG 
at the Sun position) does not change the precision on 
the reconstruction of D, but the precision on the source 
position reconstruction is better for low B,-og values. 
If we simulate the experimental resolution on energies 



by taking values such as AE/E ^ 20% instead of 
AE/E = (theoretical energies), no noticeable change 
can be seen in our results. 

V. Conclusion 

We have proposed and studied here an alternative 
way to detect individual bright UHECR sources, if 
the UHECR deflections due to the turbulent GMF are 
not too large compared to the deflections due to the 
regular component. We generated two sets of events: 
Events from a source with a l/E"^-"^ injection spectrum 
that we propagated from the source to the Earth, and 
events from some background generated according to 
the HiRes spectrum and the Telescope Array exposure. 
We mixed both sets of events, tried to detect the source 
and checked how often we got confused by detecting 
some background. 

We showed that one should look for doublets of events 
with E > 10^^'* eV, when using the depicted method. 
We studied in section |III] the dependence on several 
parameters: Unknown parameters from the magnetic 
fields, the source, the sky anisotropy at high energies and 
the experimental statistics. We showed in section HV] that 
with the parameters and assumptions adopted here, the 
precision on the source reconstruction would be domi- 
nated by the current experimental angular resolution in 
68% of cases. The local deflection power can be known 
up to 25% in 68% of cases. 

In the future, this method can be applied to experimental 
data so as to find bright UHECR proton or light nuclei 
sources. 
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